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In analogy to III-V compounds, which have significantly broadened the scope of group IV semi¬ 
conductors, we propose IV-VI compounds as isoelectronic counterparts to layered group V semi¬ 
conductors. Using ab initio density functional theory, we study yet unrealized structural phases of 
silicon mono-sulfide (SiS). We find the black-phosphorus-like ct-SiS to be almost equally stable as the 
blue-phosphorus-like /3-SiS. Both a-SiS and /3-SiS monolayers display a significant, indirect band gap 
that depends sensitively on the in-layer strain. Unlike 2D semiconductors of group V elements with 
the corresponding nonplanar structure, different SiS allotropes show a strong polarization either 
within or normal to the layers. We find that SiS may form both lateral and vertical heterostructures 
with phosphorene at a very small energy penalty, offering an unprecedented tunability in structural 
and electronic properties of SiS-P compounds. 


INTRODUCTION 

2D semiconductors of group V elements, including 
phosphorene and arsenene, have been rapidly attract¬ 
ing interest due to their significant fundamental band 
gap, large density of states near the Fermi level, and 
high and anisotropic carrier mobility [IH3]. Combina¬ 
tion of these properties places these systems very favor¬ 
ably in the group of contenders for 2D electronics ap¬ 
plications beyond graphene[5j 16] and transition metal 
dichalcogenides[7]. Keeping in mind that the scope of 
group IV semiconductors such as Si has been broad¬ 
ened significantly by introducing isoelectronic III-V com¬ 
pounds, it is intriguing to see, whether the same can be 
achieved in IV-VI compounds that are isoelectronic to 
group V elemental semiconductors. Even though this 
specific point of view has not yet received attention, there 
has been interest in specific IV-VI compounds, such as 
SnS, SnSe and GeTe, for thermoelectric and photovoltaic 
applications. [8ffT0] It appears likely that specific search 
for isoelectronic counterparts of layered semiconductors 
such as phosphorene and arsenene may guide us to yet 
unexplored 2D semiconducting IV-VI compounds that 
are stable, flexible, and display a tunable band gap. 

As a yet unexplored IV-VI compound, we study the 
layered structure of silicon mono-sulfide (SiS). We use 
ab initio density functional theory (DFT) to identify 
stable allotropes, determine their equilibrium geometry 
and electronic structure. We have identified two nearly 
equally stable allotropes, namely the black-phosphorus- 
like <a-SiS and the blue-phosphorus-like /3- SiS, and show 
their structure in Fig. [lja) and (d). Both a-SiS and f3- 
SiS monolayers display a significant, indirect band gap 
that depends sensitively on the in-layer strain. Unlike 
2D semiconductors of group V elements with the cor¬ 
responding nonplanar structure, different SiS allotropes 
show a strong polarization either within or normal to the 
layers. We find that SiS may form both lateral and ver¬ 
tical heterostructures with phosphorene at a very small 
energy penalty, offering an unprecedented tunability in 
structural and electronic properties of SiS-P compounds. 


RESULTS AND DISCUSSION 

Since all atoms in sp 3 layered structures of group V el¬ 
ements are threefold coordinated, the different allotropes 
can all be topologically mapped onto the honeycomb lat¬ 
tice of graphene with 2 sites per unit cell. m An easy 
way to generate IV-VI compounds that are isoelectronic 
to group V monolayers is to occupy one of these sites by 
a group IV and the other by a group VI element. In this 
way, we have generated the orthorhombic <a-SiS mono- 
layer structure, shown in Fig. [lja), from a monolayer of 
black phosphorene (or a-P). The hexagonal (3 -SiS mono- 
layer, shown in Fig. [ljd) , has been generated in the same 
way from the blue phosphorene (or /3- P) structure. 

The monolayer structures have been optimized us¬ 
ing DFT with the Perdew-Burke-Ernzerhof fPBEi [T2] 
exchange-correlation functional, as discussed in the 
Methods Section. We found that presence of two ele¬ 
ments with different local bonding preferences increases 
the thickness of the SiS monolayers when compared to 
the phosphorene counterparts. The 2D lattice of a- 
SiS is spanned by the orthogonal Bravais lattice vectors 
\ai\ = 4.76 A and |« 2 1 = 3.40 A, which are about 3% 
longer than the lattice vectors of black phosphorene. The 
2D hexagonal lattice of /3 -SiS is spanned by two Bravais 
lattice vectors a = \ai\ = |c? 2 | = 3.35 A, which are <i% 
longer than those of blue phosphorene. Unlike in the 
counterpart structures of the phosphorene monolayer, we 
find /3-SiS to be more stable by about 12 meV/atom than 
<a-SiS. We checked the phonon spectrum of free-standing 
<a-SiS and /3 -SiS monolayers and found no soft modes that 
would cause a spontaneous collapse of the structure. 

An apparently different layered SiS structure had been 
synthesized [13] by exposing the layered CaSi 2 structure 
to S 2 CI 2 . It appears that the S 2 CI 2 reagent reacted 
with the Ca atoms intercalated between silicene layers 
in CaSi 2 by forming CaCl 2 and saturating the valencies 
of the silicene layers by — S—S bridges. Even though no 
structural information has been provided in that study, 
we identified a likely candidate structure and found it 
to be less stable than the postulated SiS allotropes, on 
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FIG. 1. (Color online) Atomic and electronic structure of [(a) to (c)] o-SiS and [(d) to (f)] /3-SiS monolayers, [(a) and (d)] 
Ball-and-stick models of the geometry, with S and Si atoms distinguished by size and color and the Wigner-Seitz cell indicated 
by the shaded region, [(b) and (e)] The electronic band structure and [(c) and (f)] the electronic density of states (DOS) of the 
systems. The energy range between Ep and 0.2 eV below the top of the valence band, indicated by the green shading, is used 
to identify valence frontier states. 


which we focus next. 

Since the electronegativity of S is higher than that of 
Si, we expect an electron transfer from Si to S atoms. 
Based on a Mulliken population analysis, we estimate a 
net transfer of 0.3 electrons in a-SiS and 0.2 electrons in 
/3-SiS from silicon to sulfur atoms. This charge redistri¬ 
bution, combined with the nonplanarity of the structure, 
causes a net polarization. We find an in-plane polariza¬ 
tion for <a-SiS and and an out-of-plane polarization for 
/3-SiS. 

Whereas DFT generally provides an accurate descrip¬ 
tion of the total charge density and equilibrium geometry, 
interpretation of Kohn-Sham energy eigenvalues as quasi¬ 
particle energies is more problematic. Still, we present 
our DFT-PBE results for the electronic structure of a- 
SiS and /3-SiS monolayers in Fig. [l] Even though the fun¬ 
damental band gaps are typically underestimated in this 
approach, the prediction that both systems are indirect- 
gap semiconductors is likely correct. Our calculated band 
structure and the corresponding density of states for a- 


SiS, presented in Fig.[IJb,c), suggest that the fundamen¬ 
tal band gap value E g = 1.44 eV should be significantly 
larger than in the isoelectronic black phosphorene coun¬ 
terpart. The band structure near the top of the valence 
band shows a significant anisotropy when comparing the 
T — X and T — Y directions. In analogy to phospho¬ 
rene, <a-SiS should exhibit a higher hole mobility along 
the x-direction than along the ^/-direction. 

The DFT-based fundamental band gap E g = 2.26 eV 
in the monolayer of /3-SiS is even larger. The band 
structure in the symmetric honeycomb lattice of /3-SiS 
is rather isotropic, as seen in Fig. [lje). The top of the 
valence band is very flat, resulting in a heavy hole mass 
and a large density of states (DOS) in that region, as 
seen in Fig. []Jf). 

The character of frontier states is of interest not 
only for a microscopic understanding of the conduction 
channels, but is also crucial for the design of optimum 
contacts. m Whereas DFT-based band gaps are typi¬ 
cally underestimated as mentioned above, the electronic 
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FIG. 2. (Color online) Electron density p v b associated with 
states in the energy range between the Fermi level Ef and 
0.2 eV below the top of the valence band, shown by the shaded 
region in Fig. [I] in a monolayer of (a) o-SiS and (b) /3-SiS. 
p v b — 2.5xlO - ' T e/A 3 contours are superposed with ball-and- 
stick models of related structures. 


structure of the valence and the conduction band region 
in DFT is believed to closely correspond to experimen¬ 
tal results. In Fig. [2| we show the charge density associ¬ 
ated with frontier states near the top of the valence band. 
These states, which correspond to the energy range high¬ 
lighted by the green shading in the band structure of 
<a-SiS in Fig. [ljb) and that of /3-SiS in Fig. [lje), cover 
the energy range between the mid-gap and 0.2 eV be¬ 
low the top of the valence band. The valence frontier 
states of <a-SiS in Fig. §a) and (3 -SiS in Fig. |2jb) are 
similar in spite of the notable charge density differences 
caused by the larger DOS of ft -SiS in this energy range. 
The difference between S and Si atoms is also reflected in 
the character of the valence frontier states at these sites. 
Whereas p z states contribute most at sulfur sites, silicon 
sites contribute a mixture of s and p z states. These fron¬ 
tier states differ from those of phosphorene, which are 
related to lone pair electron states. 

Similar to phosphorene, the fundamental band gap 
value of SiS also depends sensitively on the in-layer 
strain, as seen in Fig. [3j Due to their nonplanarity, 
accordion-like in-layer stretching or compression of SiS 
structures may be achieved at little energy cost as shown 
in the Supporting Information. The energy cost is partic¬ 
ularly low for a deformation along the soft ^-direction, re¬ 
quiring <20 meV/atom to induce a ±9% in-layer strain. 
We believe that in view of the softness of the structure, 
similar strain values may be achieved during epitaxial 
growth on particular incommensurate substrates. We 
also note that tensile strain values such as these have been 
achieved experimentally in suspended graphene mem¬ 
branes that are much more resilient to stretching due to 
their planar geometry and stronger bonds fl5l[l7] . Con¬ 


sequently, we believe that strain engineering is a viable 
way to effectively tune the fundamental band gap in these 
systems. 

Our results for <a-SiS in Fig.[3ja) indicate that the band 
gap decreases when the structure is compressed and in¬ 
creases slightly when it is stretched. The largest change 
in the band gap, namely its reduction to 0.3 eV, may be 
achieved during a 9% compression. As seen in Fig. |3|b), 
we expect the fundamental band gap of /3-SiS to be re¬ 
duced during both stretching and compression. Within 
the ±9% range, we find that the band gap may be tuned 
in the range from ^1.0 —2.3 eV. This high degree of band 
gap tunability in SiS appears particularly attractive for 
potential applications in flexible electronics. 

Since the geometry and lattice constants of SiS and 
phosphorene are very similar, it is likely that the two 
could interface naturally in lateral and vertical het¬ 
erostructures, thus further advancing the tunability of 
their electronic properties. In Fig [4j we present elec¬ 
tronic structure results for bilayers consisting of an SiS 
and a phosphorene monolayer in both a- and /3-phases 
as the simplest examples of vertical heterostructures. We 
have optimized the bilayer structures assuming commen- 
surability, i.e. setting the primitive unit cells of each 
monolayer to be the same. The optimum geometry of 
the <a-(SiS/P) bilayer is shown in Fig Qa) and that of 
the /3- (SiS/P) bilayer in Fig[4jb). 

We find the inter-layer interaction in the two bilayer 
systems to be rather weak, amounting to <20 meV/atom 
based on our DFT-PBE calculations. Taking better ac¬ 
count of van der Waals interactions would likely increase 
the calculated interaction energy and reduce the large 
interlayer separation of d<5.4 A. Whereas the precise 
interlayer interaction and separation are not of primary 
concern here, our most important finding is that the weak 
interaction is not purely dispersive in nature, since we 
find a substantial rehybridization of states between the 
adjacent SiS and phosphorene layers. Consequently, as 
we show more explicitly in the Supporting Information, 
the bilayer band structure is not a mere superposition of 
the two monolayer band structures in the same assumed 
geometry. 

As shown in Fig.[4^b), the <a-(SiS/P) bilayer is a direct- 
gap semiconductor with £^0.8 eV, smaller than the 
band gap in either isolated monolayer. Since the gap 
is indirect in the <a-SiS monolayer, whereas it is direct in 
<a-P, the cause for the direct band gap in <a-SiS/P may 
be the dominance of P states near the gap. We find this 
conjecture confirmed by examining the frontier orbitals 
of the bilayer. As we show in the Supporting Informa¬ 
tion, the frontier states in the conduction band region of 
the bilayer are essentially purely P-based, whereas those 
in the valence band region have only a small contribution 
from the SiS layer. 

The electronic band structure of the AA-stacked /3- 
(SiS/P) bilayer is presented in Fig. |4jd). The funda- 
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FIG. 3. (Color online) Electronic band gap of (a) o-SiS and (b) /3-SiS monolayers as a function of the in-layer strain. Besides 
the uniform strain, we also show results for uniaxial strain along the x and y direction in the anisotropic structure of o-SiS in 
(a). The lines are guides to the eye. 

(a) a-(SiS/P) bilayer (c) p-(SiS/P) bilayer 



FIG. 4. (Color online) Optimum geometry and electronic band structure of [(a) and (b)] an a-(SiS/P) bilayer and [(c) and (d)] 
a f3-( SiS/P) bilayer. The optimum stacking of the SiS and the phosphorene monolayers in the a-(SiS/P) bilayer in (a) is AB 
and that in the (3- (SiS/P) bilayer is AA. 


mental band gap is indirect in this system and its value 
Egtt 1.7 eV is again smaller than in isolated blue phos¬ 
phorene and (3 -SiS monolayers. As seen in the Support¬ 
ing Information, also the frontier states of (3- (SiS/P) are 
dominated by P orbitals both in the valence and the con¬ 
duction band region. 

Since both SiS and phosphorene are rather flexible, 
they may adjust to each other and form also in-layer 
heterostructures at little or no energy penalty. We con¬ 


structed two types of SiSP 2 heterostructures for both 
a and (3 allotropes and show their geometry and elec¬ 
tronic structure in Fig.|5j One type of heterostructures, 
shown in Fig. [5ja) and[5fe), contains P-P and Si-S atom 
pairs completely separated from like atom pairs. The 
other type of heterostructures, presented in Fig. [5jc) 
and[5][g), contains alternating, contiguous SiS and phos¬ 
phorus chains. These structures maintain a rectangu¬ 
lar lattice with four atoms per unit cell. Generally, we 
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FIG. 5. (Color online) Optimum geometry and electronic band structure of [(a) to (d)] two different in-layer heterostructures 
of <a-SiSP 2 and [(e) to (h)] two different in-layer heterostructures of /3-SiSP2- The heterostructures differ in the arrangement 
of P atoms. 


found the in-layer heterostructures to be less stable than 
pure phosphorene and SiS monolayers. The least sta¬ 
ble heterostructures among these are those with isolated 
P-P or Si-S atom pairs, shown in Fig. |5|a) and [5je) , 
which are ~0.2 eV/atom less stable than SiS and phos¬ 
phorene monolayers due to their highly frustrated ge¬ 
ometries. The heterostructures with contiguous phos¬ 
phorus and SiS chains, presented in Fig. §c) and |5jg), 
may better optimize the nearest-neighbor environment. 
This causes less frustration, making these systems only 
~0.1 eV/atom less stable than isolated SiS and phospho¬ 
rene monolayers. 

We found all four in-layer heterostructures to be 
indirect-gap semiconductors. As in the vertical het¬ 
erostructures, we found the fundamental band gaps to 
be substantially smaller than in isolated SiS and phos¬ 
phorene monolayers. As seen in Fig. |5][a,b) and §e,f), 
the fundamental band gap E g is close to 0.5 eV in the 
less stable heterostructures with isolated Si-S and P-P 
pairs. We found larger band gap values in the more sta¬ 
ble heterostructures with contiguous SiS and P chains, 
namely E g & 0.6 eV in <a-SiSP 2 shown in Fig. [5^c,d) and 
E g tt 1.1 eV in f3- SiSP 2 shown in Fig. [5jg,h) , both nearly 


half the value of the smaller band gap in either SiS or 
phosphorene monolayers. These findings indicate an in¬ 
triguing possibility of isoelectronic doping as an effective 
way to tune the electronic properties of SiSP n systems. 

In view of the different successful synthesis approaches 
used to form layered IV-VI compounds such as sis mu or 
SnS[18j with a structure similar to <a-SiS, we believe that 
a suitable synthesis path will be found to form also the 
allotropes introduced in this study. The weak inter-layer 
interaction in layered pure SiS compounds or SiS/P het¬ 
erostructures should allow for a mechanical exfoliation 
yielding monolayers and bilayers, same as in graphene 
and phosphorene. Chemical Vapor Deposition (CVD), 
which had been used successfully in the past to grow 
graphene [T9l [20] and silicene[2T], may become ultimately 
the most common approach to grow few-layer SiS on spe¬ 
cific substrates. 

As suggested above, compounds that are isoelectronic 
to group V layered semiconductors are not limited to pure 
IV-VI systems, but may contain other group V elements 
in the same layer or in a vertical few-layer heterostruc¬ 
ture. The initial interest in black phosphorene could thus 
be significantly expanded by considering the whole range 































6 


of group IV elements including Si, Ge, Sn, and Pb, group 
VI elements such as S, Se, and Te, and group V elements 
such as P, As, Sb and Bi, leading to a virtually limitless 
number of compounds and structural phases. The large 
family of IV-VI systems that are isoelectronic to group 
V elemental semiconductors will include 2D semiconduc¬ 
tors with a sizeable fundamental band gap, a high carrier 
mobility and chemical stability. Expanding this family to 
heterostructures containing also group V elements should 
provide a way to effectively tune the electronic properties 
of the pristine structures, which will likely bring unlim¬ 
ited richness to the field of 2D semiconductors. 

In conclusion, we have proposed IV-VI compounds 
as isoelectronic counterparts to layered group V semi¬ 
conductors in analogy to III-V compounds, which have 
significantly broadened the scope of group IV semicon¬ 
ductors. Using ab initio density functional theory, we 
have identified yet unrealized structural phases of silicon 
mono-sulfide (SiS) including the black-phosphorus-like a- 
SiS and the almost equally stable, blue-phosphorus-like 
/3- SiS. We found that both a-SiS and /3 -SiS monolayers 
display a significant, indirect band gap that depends sen¬ 
sitively on the in-layer strain. Unlike 2D semiconductors 
of group V elements with the corresponding nonplanar 
structure, different SiS allotropes are polar. We found 
that SiS may form both lateral and vertical heterostruc¬ 
tures with phosphorene at a very small energy penalty, 
offering an unprecedented tunability in structural and 
electronic properties of SiS-P compounds. Combining 
other group IV and group VI elements with group V el¬ 
ements is expected to lead to a large family of layered 
semiconductors with an unprecedented richness in struc¬ 
tural and electronic properties. 

METHODS 

Our computational approach to study the equilibrium 
structure, stability and electronic properties of SiS struc¬ 
tures is based on ab initio density functional theory 
with the Perdew-Burke-Ernzerhof (PBE)[12] exchange- 
correlation functional as implemented in the SIESTA [ 2 2] 
and VASP [23H26j codes. We use periodic boundary 
conditions throughout the study, with few-layer struc¬ 
tures represented by a periodic array of slabs sepa¬ 
rated by a vacuum region of >15 A. In our SIESTA 
calculations we used norm-conserving Troullier-Martins 
pseudopotentials[27] and a double-^ basis including po¬ 
larization orbitals. The reciprocal space has been sam¬ 
pled by a fine grid [28] of 8x8x1 ^-points in the Brillouin 
zone of the primitive unit cell of the 2D structures. We 
used a mesh cutoff energy of 180 Ry to determine the 
self-consistent charge density, which provided us with a 
precision in total energy of <2 meV/atom. In VASP 
calculations, we used an energy cut-off of 500 eV and 
the same k-point sampling as mentioned above. All ge¬ 


ometries have been optimized using the conjugate gra¬ 
dient method [29], until none of the residual Hellmann- 
Feynman forces exceeded 10“ 2 eV/A. Equilibrium struc¬ 
tures and energies based on SIESTA have been checked 
against values based on the VASP code. 
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